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Abstract 

O ' 

CN ■ We derive bounds on the mixing between the left-chiral ("active") and the right- 

ed , chiral ("sterile") neutrinos, provided from the combination of neutrino oscillation data 

and direct experimental searches for sterile neutrinos. We demonstrate that the mixing 
of sterile neutrinos with any flavour can be significantly suppressed, provided that the 
angle #13 is non-zero. This means that the lower bounds on sterile neutrino lifetime, 
coming from the negative results of direct experimental searches can be relaxed (by as 
much as the order of magnitude at some masses). 

We also demonstrate that the results of the negative searches of sterile neutrinos with 
PS191 and CHARM experiments are not applicable directly to the see-saw models. The 
reinterpretation of these results provides up to the order of magnitude stronger bounds 
on sterile neutrino lifetime than previously discussed in the literature. We discuss the 
implications of our results for the Neutrino Minimal Standard Model (the z^MSM). 
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1 Introduction 



Transitions between neutrinos of different flavours (see e.g. [l] for a review and Refs. 0, 0] for 
the recent update of experimental values) are among the few firmly established phenomena 
beyond the Standard Model of elementary particles. The simplest explanation is provided by 
the "neutrino flavour oscillations" - non-diagonal matrix of neutrino propagation eigenstates 
in the weak charge basis. While the absolute scale of neutrino masses is not established, 
particle physics measurements put the sum of their masses below 2 eV [i| while from the 
cosmological data one can infer an upper bound of 0.58 eV at 95% CL j^]. 

A traditional explanation of the smallness of neutrino masses is provided by the see-saw 
mechanism It assumes the existence of several right-handed neutrinos, coupled to their 

Standard Model (SM) counterparts via the Yukawa interaction, providing the Dirac masses, 
M D , for neutrinos. The Yukawa interaction terms dictate the SM charges of the right-handed 
particles: they turn out to carry no electric, weak and strong charges; therefore they are often 
termed "singlet," or "sterile" fermions. Sterile neutrinos can thus have Majorana masses, M s , 
consistent with the gauge symmetries of the Standard Model. If the Majorana masses are 
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much larger than the Dirac ones, the type I seesaw formula [6-9] holds, expressing the mass 
matrix of observed neutrinos (Ai u ) via 

M v = M D M S 'Ml (1) 

where M. u is a 3 x 3 matrix of active neutrino masses, mixings, and (possible) CP- violating 
phases. The masses of sterile neutrinos are given by the eigenvalues of their Majorana mass 
matrix (with the corrections of the order M^/Mg). They are much heavier than the active 
neutrino masses as a consequence of ([T]). 

Numerous searches for sterile neutrinos in the mass range up to ~ 100 GeV had been 



performed (see the corresponding section in Particle Data Group |4([_!] see also [10J and refs. 
therein). These searches provided upper bounds on the strength of interaction of these neutral 
leptons with the SM neutrinos of different flavours - active-sterile neutrino mixing angles 



Md,, 



for sterile neutrino with the mass M s % These bounds then can be interpreted as 



lower bounds on the lifetime of sterile neutrinos r, via 



_ 1= GlM[ 
Ts " 96vr 3 



£C4 a) , (2) 



where the sum runs over various particles to which sterile neutrino can decay, depending on 
their mass (u, e ,fx ,r , tt, K, heavier mesons and baryons) and dimensionless quantities B\ 
depend on the branching ratios (see Appendix|X]for details). The lower bound on the lifetime 
t s is usually dominated by the least constrained mixing angle, ^ (as will be shown later). 

This bound can be made stronger if one assumes that the same particles are also responsible 
for the neutrino oscillations. The see-saw formula ([I]) allows to fix (at least partially) the ratio 
of the mixing angles i^/^i of different flavours. In the simplest case when only two sterile 
neutrinos are present (the minimal number, required to explain two observed neutrino mass 
differences) the ratio of all three mixing angles can be fixed, which makes the bounds on the 
lifetime significantly stronger (essentially, they are determined now by the strongest, rather 
than the weakest direct bound on •&„). When confronted with the cosmological bounds (e.g. 
from Big Bang Nucleosynthesis [O, they seem to close the window of parameters for 



sterile neutrinos with the mass lighter than about 150 MeV [131 . [14 . 

In this paper we demonstrate that the bounds on sterile neutrino lifetime from accelerator 
searches can be relaxed even in the case when there are only two sterile neutrinos, responsible 
for the observed neutrino oscillations. Namely, we show that the recent neutrino oscillation 
data [3] allow for such a choice of the active-sterile Yukawa couplings that the mixing of sterile 
neutrinos with any given flavour can be strongly suppressed (provided that the active neutrino 



mixing angle (9 13 0, as hinted recently by MINOS 15] and T2K [16| experiments). This 
changes the restrictions on the sterile neutrino lifetime in some cases by as much as 10 times 
in the normal hierarchy of active neutrino masses and two times in the inverted hierarchy. The 
results of this paper partially overlap with [13 (also 0), and we compare in the corresponding 
places to the previous works. 



J littp : //pdglive . lbl . gov/Rsummary .brl?nodein=S077&inscript=Y 

2 Here and below we use the letter i? for active-sterile mixing angles (denned by Eq. (1111) below) while 
reserving #12, #13 and #23 for the measured parameters of the active neutrinos matrix M. v . These quantities 
•d a are often denoted |T4a| 2 or It^al 2 in the experimental papers, to which we refer. Here and below the Greek 
letters a, /3 are flavour index e, [A, r and i,j = 1, 2, 3 denote active neutrino mass eigenstates. 
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The paper is organized as follows: in Section [2] we briefly describe the type I see-saw model 
we use. We then analyze the see-saw equations and demonstrate that the mixing with any 
flavour d 2 a can become suppressed and investigate relations between different mixing angles 
imposed by the see-saw mechanism (Section [3]). Section H] is devoted to the overview of the 
experiments, searching for sterile neutrinos with the masses below 2 GeV, and the way one 
should interpret their results to apply to the see-saw models that we study. Section summa- 
rizes our revised bounds on mixing angles and translates them into the resulting constraints on 
sterile neutrino lifetime (Figs. [6]). We conclude in Section El discussing possible implications 
of our results. 



2 Sterile neutrino Lagrangian 

The minimal way to add sterile neutrinos to the Standard Model is provided by the Type I 
see-saw model p-0] (see also [l9[ and refs. therein): 

AC N = iNjd^Nj - (F aI L a N& ~ NfNj + h.A , (3) 

i,j=i \ - / 

where F a i are new Yukawa couplings, $j is the SM Higgs doublet, $j = This model 

is renormalizable, has the same gauge symmetries as the Standard Model, and contains Af 
additional Weyl fermions Nj - - sterile neutrinos (Nf being the charge-conjugate fermion, 
in the chiral representation of Dirac 7- matrices Nf = i^ 2 N\). The number of these singlet 
fermions must be M > 2 to explain the data on neutrino oscillations. In the case of M = 2 
there are 11 new parameters in the Lagrangian (J3j), while the neutrino mass/mixing matrix 
M. v has 7 parameters in this case. The situation is even more relaxed for M > 2. The 
see-saw formula (pQ) does not allow to fix the scale of Majorana and Dirac M D Q j- = F a j(Q) 
masses. In this work we will mostly concentrate on sterile neutrinos with their masses M s in 
the MeV-GeV range - the range that is probed by the direct searches. To further simplify 
our analysis we will concentrate on the case when the masses of both sterile neutrinos are 
close to each other (so that AM <C M s ). One important example of such model is provided 



by the Neutrino Extended Standard Model (the z/MSM) ([2fJ, |21|, see [14| for review). The 
z/MSM model contains 3 sterile neutrinos, whose masses are roughly of the order of those 
of other leptons in the Standard Model. Two of these particles (approximately degenerate 
in their mass) are responsible for baryogenesis and neutrino oscillations and the third one is 
playing the role of dark matter. The requirement of dark matter stability on cosmological 
timescales makes its coupling with the Standard Model species so feeble, that it does not 
contribute significantly to the neutrino oscillation pattern |22j. Therefore, when analyzing 
neutrino oscillations, the Ni can be omitted from the Lagrangian and index / in the sums 
runs through 2 and 3 only: 

£ see . saw = C SM + Niidp-fNj - M D , aI u a Nj - M* j(xI N lVa - M S N<N 3 - M S N 3 N< . (4) 



This parametrization coincides with 23j, Eq. (2.1)]. Note that we use basis of singlet neutrinos 



where the mass matrix is off-diagonal. Recent computation of the baryon asymmetry of 



the Universe in the z/MSM 241 demonstrated that sterile neutrinos with the mass as low as 



several MeV can be responsible for baryogenesis and neutrino oscillations within the z/MSM. 
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This prompts us to re-analyze the implication of negative direct experimental searches for 
the Yukawa couplings of sterile neutrinos with the MeV masses. We limit our analysis by 
M s < 2 GeV, as for the higher masses the existing experimental bounds do not probe the 
region of mixing angles, required to produce successful baryogenesis in the z/MSM 24 . 



3 Solution of the see-saw equations 

In this Section we investigate how mixing angles between active and sterile neutrinos are 
related to parameters of the observable neutrino matrix M. v . We will demonstrate that the 
mixing angle d 2 in the case of normal hierarchy and the mixing angles "d 2 or d 2 in the case of 
inverted hierarchy, can become suppressed as we vary the parameters of the neutrino matrix 
away from their best-fit values (but within the experimentally allowed 3a bounds). 



3.1 Parametrization of the Dirac mass matrix 

We use the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) parametrization of the neutrino ma- 
trix M v (see e.g. Eqs.(2.10) and (2.12) of fl) 

M v = V* diag(m ie - 2 S m 2 e~ 2i ^ m 3 )V f , (5) 

where V is the unitary matrix, whose explicit form is reminded in Appendix [Bj Redefining a 
Dirac mass matrix ad_| 

M D — y M D = V T M D , (6) 
we can rewrite the see-saw relation ([1]) in the following form: 

diagKc-* m 2 e- 2 *, m») y = . ^^+^^ , ( 7 ) 

The rank of the active neutrino mass matrix A4 U is 2 in the case of two sterile neutrinos, 
meaning that one of the masses rrii is zero. Two choices of "hierarchies" of the mass eigenvalues 
are possible. The first one is called normal hierarchy (NH) and corresponds to < m\ < m 2 < 
m 3 . The second one is called inverted hierarchy (IH) and is realized for < m 3 < mi < m 2 . 

Once the mass M s is fixed, the solutions of Eq. (J7|) contain one unknown complex pa- 
rameter, z. Its presence reflects a symmetry of the see-saw relation (J7j) under the change 
(M D ,i2, M D! is) — > (zM D) i2, z" 1 M-D,a) [25(- It is this freedom that does not allow to fix the 
absolute scale of M D (i.e. the value of d 2 ) even if M s is chosen. 

The change z — > z~ l is equivalent to the redefinition of N 2 — >■ N 3 , N 3 — > N 2 together with 
shift of the Majorana phase £ — > £ + 7r in ([7]) . Therefore in subsequent analysis we will choose 
\z\ > 1 without the loss of generality. 

3.2 Normal hierarchy 

For normal hierarchy the explicit see-saw relation is 

diag(0, m 2 e-*«, ,n 3hl = + . (8) 



3 The Dirac matrix M D has indexes / = 2, 3 and i,j = 1, 2, 3 - neutrino propagation basis 
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Normal hierarchy 



Inverted hierarchy 



Am 2 21 


(7.09 


-8.19) x 10" 


" 5 eV 2 








Am 2 31 


(2.14 


- 2.76) x 10- 


" 3 eV 2 


Am 2 3 


(2.13 


- 2.67) x 10~ 3 eV 2 


sin 2 9 12 


0.27- 


-0.36 










sin 2 6> 23 


0.39- 


-0.64 










sin 2 6»i3 


0.001 


-0.035 




sin 2 6»i3 


0.001 


- 0.039 



Table 1: The 3a bounds on the parameters of the mass matrix M. v , adopted from |3|. Here 



Am. 



mi — m- 



The boundaries for inverted hierarchy as the same as for the normal one, 



unless written explicitly. 

Diagonal components of this matrix equation give 

M d ,i 2 Md,i3 = 0, M D , 22 M D , 23 = X -m 2 M s t- 2i ^ M D , 32 M D ,33 = \m z M a . (9) 

Using m 2 ,m3 ^ we find that M Di22 , M D>2 3, M Dj 3 2 , M Dj 33 are all non-zero. Analysis of non- 
diagonal terms reveals that both M D ,i2 an d Md,13 are zero and there are two general solutions 
(c.f. H): 



M^2 



iz 



M, 



(0, ±ie 



™3j, 



iz 



(0, =p ie-^V^, v 7 ^) • (10) 



The solution M+ with £ = ip + ir equals to M~ with £ = ip. It allows us to consider only one 
solution on the interval < £ < 2tc. In what follows we therefore omit the superscript +Q 
The mixing angles of the active-sterile neutrinos are defined as follows: 



2tf 



E KMoiicvr = E ic^aUcvi 2 = E i^w 



Inserting the explicit solution ffTOl) for M D results in 



M 2 



(11) 



4M S 



AM s \z\ 2 



(12) 



For \z\ S> 1 the contribution of M D>i3 is suppressed compared with that of M D] j 2 and therefore 
we neglect the former (we will comment below on the case \z\ > 1). 

As the value of the Majorana phase £ is undetermined experimentally, the condition d a = 
is satisfied iff m^V^I 2 = m 2 |V a2 | 2 (we neglect second term on the r.h.s. of ffT2]) ). For the 
electron flavour (a = e) it translates into 



sin 6*12 — = tan #13, 
m 3 



(13) 



4 Unlike the parametrizations used e.g. in Ref. 17, 2|| 26 1 this way of parametrizing the solution of the 
see-saw equations shows that there is only one branch of solutions, with all other related to it via redefinitions 
AT 2 -H> A*3 and shift of the Majorana phases. In particular in the parametrization we used it is much easier to 
analyze whether mixing angles become zero. The relation \z\ — exp(Imw) holds, where the parameter lu was 
employed in [17j . 
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which, in principle, can be satisfied only for non-zero #13. This result has been already obtained 
in 



17 



The bounds on the parameters of the mass matrix Ai u at the 3a level that we use are 
shown in Table [TJ Note that in this paper we do not take into account statistical correlations 
between different oscillation parameters, allowing them to vary independently within their 
3a intervals. Consequently, we obtain the 3a intervals for the combinations of parameters, 
entering Eq. (ITBl : 

0.043 < sin 2 fl 12 ^ < 0.070, 

0.001 < tan 2 fl 13 < 0.035. 1 ; 

It means that the relation (fl3|) does not hold for the neutrino oscillation parameters, presented 
in Table [TJ Therefore the mixing angle °d\ cannot become zero, but has a non-trivial lower 
bound. To find the minimal value that it can reach, we consider the ratio of the angles 
$ 2 / {Q\ + $ 2 + i? 2 .). Due to the unitarity of V, the denominator is 

w 2W2 £ V:,M^ 2 V ai Ni; a2 = _-L^|M Di/32 | 2 = + m 3 ). (15) 

a s a,pn s /3 s 

Let us denote the ratio of the mixing of sterile neutrinos with one flavour to the sum of all 
mixings by T a , 

T a = (16) 



Then we get the following expression for T e : 



Te = ^ • (17) 

The minimum is achieved if we push #12 and Am^ to their 3a lower boundaries, #13 and 
Am^ to their upper boundaries, and choose £ = —71/ 2. The maximum is achieved when we 
set Amjj equal to its lower bound, #13, 612 and Am^i to their upper bounds, and by choosing 
the Majorana phase a = tt/2. The bounds on T e from Table [2] translate into the bounds for 
the muon and tau flavours combined: 

0.83 < T u + T T < 0.9997. (18) 



The minimum and maximum of different T a are listed in the Table [2]1_| This analysis was 
conducted in approximation of large \z\. See Appendix [C] for the account of finite- 1 z \ effects. 

3.3 Inverted hierarchy 

Similarly to the previous case, for the inverted hierarchy we get a solution of the see-saw 
equations (J7J) 



M D)i2 = iz\l^(e~^^rT l , ie-*V™* °)> ^3 = J ^(e^V™^, -ie*y/m, 0) (19) 



5 Notice that if one uses the results of [27} instead of @, Q to determine the upper and lower bounds on the 
neutrino mass matrix, the results, listed in the Table[3]can get modified (mostly, because the Ref. [2lJ provides 
broader interval for the possible values of 6*13). For example, it becomes possible to suppress completely rJ e fl7j j . 
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Normal hierarchy 


Inverted hierarchy 


1.7 x 10~ 3 < T e < 0.15 
0.09 < Tjj < 0.89 
0.08 < T r < 0.88 


0.02 <T e < 0.98 
3 x 10~ 5 <T^< 0.60 
1.5 x 10~ 5 <T T < 0.62 


The ranges are based on 2a bounds 



Normal hierarchy 


Inverted hierarchy 


3 x KT 4 < T e < 0.17 
0.07 < T M < 0.92 
0.06 <T T < 0.90 


0.02 < T e < 0.98 
< T M < 0.63 
< T T < 0.65 


The ranges are based on 3a bounds 



Table 2: The ratio of the sterile neutrino mixing with a given flavour a to the sum of the 
three mixings, T a (defined by (TIB"]) ). Left table shows the upper and lower values of T a when 
parameters of neutrino oscillations are allowed to vary within their 2a boundaries (taken from 
[3]). The right table shows the results when the parameters of active neutrino oscillations are 
varied within their 3a limits (see Tabled]). See the footnote [5] regarding the use of different 
3a bounds. 



for < £ < 2n. In this case ft 2 or -d 2 . can become very suppressed, as we will show soon. 
The mixing angles are 

^ = 4M I ^™ lVal ~ ie ^~° V^V a2 1 2 + j^j-r 2 | VW«1 + ie l{ ^° s/m^V a2 1 2 . (20) 



For \z\ ^> 1 they can become close to zero only if y^rni|T^:i| = v^m^lV^I- For a = /x this 
condition translates into 

Itan^a + sin^atan^ae^l = ./— 11 - sin0i 3 tan0 12 tan0 23 e~^|. (21) 

V mi 

For the parameter set close to the best fit, left-hand side is less than the right-hand side, 
because then sin# 13 w 0, while tan^ 12 < 1 and mi ~ m 2 . To attain the equality one has to 
push left-hand side up and the right-hand side down. = makes phases of both complex 
terms inside |...| on the left-hand side equal, thereby the absolute value of their sum becomes 
maximal. Simultaneously the right-hand side becomes minimal. For this specific choice of the 
Dirac angle the equality ( 12 1]) turns into 

,/^-tan# 12 

\ mi L * 

v - sin0 13 tan0 23 . (22) 



^tan0 12 + l 

The 3a bounds for inverted hierarchy in general are the same as for the normal one (see 
Table [1]) with the exception of the "atmospheric" mass difference and #i 3 bound, that slightly 
differ. Using these values we find 

tan flu 

0.14 <-^= < 0.24, 0.03 < sin 0i 3 tan 23 < 0.26. (23) 

,M tan 0i2 + 1 

We see that two regions overlap, therefore the relation (j2"Tl) can be satisfied and ft 2 can be zero 
in a wide region of values of the parameters of the neutrino oscillation matrix. See, however, 
Sec. 13.41 below. 
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Similarly, the condition d T = (for <p = tt) translates into 



- tan U 

mi " 



'^tan0 12 + 1 

mi lz 



sin 0i3 cot 023, (24) 



and can be satisfied, because the quantity on the right hand side varies from 0.02 to 0.25, well 
within the range of d23])@ 

On the other hand, $ e can be zero only if 

cot 0i 2 = ,/— (25) 
V mi 

can be realized. The left hand side is always larger than the right hand side (within the 3cr 
region), therefore no d e suppression can occur. However it is important to know what minimal 
value this mixing angle can reach. According to Eq. (l20|) electron mixing angle is given by 

U|2 

T?e = -TJ7 C0S2 C ° S2 ^ 12 + 1712 Sin2 ^ 12 + Sil1 ^ _ ^ Sil1 2 ^12A/ m l m 2) • (26) 

For C, — C = — 7r/2 this quantity is minimal 



^ e ,mm = c os 0i3 ( v / ^i"cos0i 2 - y/m^sm8 12 ) . (27) 
To compare it with the other mixing angles, we note that the relation 

12 



AM, 

a 

holds (similar to Eq. fll5j) in the case of normal hierarchy). Therefore 



®l,min COS 2 013 ( fmi . V . . 

COS 012 — \ SH1012 . (29) 



EC l + ^V Vmi 



The results of the analysis are listed in Table [2J From the upper bound on T a we derive 
the bound 

T M + T T > 0.02 . (30) 

We see that in this mass hierarchy it is possible for the overall coupling of the sterile neutrino 
to both fi and r flavours to become tiny compared to the electron flavour coupling. 



6 It was pointed out in [17( that both z9 M and $ T can be suppress in inverted hierarchy, for #13 = 0. For 

this to happen the relation y^j^ = tan 6*12 should hold, as one can also see from Eqs. and AMI) . The 

corresponding value of #12 is however well outside the 3<r interval. The general case #13 ^ has not been 
analyzed in [17| . 
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Flavour a 


(^) min @ 1 MeV 


z 


e 


10" 9 


1.7 




1(T 8 


1.5 


(a) NH, best-fit 


Flavour a 


(^) min @ 1 MeV 


\z\ 


e 


io- 8 


2.3 




2 x 10~ 9 


3.5 



Flavour a 


(^) min @ 1 MeV 


\z\ 


e 


2 x 1(T 10 


5 




6 x 1CT 9 


1.8 


(b) NH, 3a 


Flavour a 


(^) min @ 1 MeV 


\z\ 


e 


6 x 1(T 9 


2.7 




see text 





(c) IH, best-fit 



(d) IH, 3a 



Table 3: Minimal values of the active-sterile mixing angles i?^, obtained using the best fit 
values of neutrino oscillation parameters (left columns) or by varying the neutrino oscillation 
data within their 3a intervals, listed in Tabled] (right columns). The values for (f^mm are 
provided for sterile neutrinos with the mass M s = 1 MeV. For other masses one should 
multiply the listed values by ( MeV/M s ). Columns show the values of |z| for which the 
minimum in fl3TD is reached. 



3.4 Minimal mixing angles in the z/MSM 

Finally, we find the minimal values of the sterile neutrino mixing angles in the z/MSM, com- 
patible with the neutrino oscillation data. These angles will turn out to be much smaller than 
the experimental upper bounds in all regions of masses, probed by the experiments. A general 
solution of the see-saw equations (|T2l . fT20|) gives ^ as a function of \z\: 

i}l = A a \z\ 2 + h (31) 
\ z \ 

with coefficients A a and B a independent of \z\. The minimum of this expression is reached 
for \z\ 2 a = a/ B a /A a < 1 and is given by 

(^) min = 2y/A a B a . (32) 

To find an absolute lower bound on the mixing angle for a given sterile neutrino mass, we vary 
this expression over the parameters of neutrino oscillations. The resulting mixing angles and 
the corresponding values of \z\ are listed in Table [3lPI 

For the mixing angles q d 2 flT in the case of inverted hierarchy, different consideration takes 
place. Formally for infinitely large \z\ they would become zero. The value of |z|, however, is 
bounded from above from the experiments \z\ < z u \,. To see this, one should look at Fig. [TJ 
where the resulting bounds are given. Therefore couplings to /i ant r neutrinos are finite for 
any phenomenologically plausible parameter set of the z/MSM. The appropriate estimates of 
mixing angles can be provided for B a given by L 1 ^ 1 (I46|I47I) . along with A^ T — 0, z = z ub 

<r > 2 x KT 8 v^r, (IH) (33) 



7 Notice, that the ratios of the mixing angles do not reach their minima/maxima when the angles 

are minimized. The "interesting" values of \z\ (those that can relax the bounds on the lifetime) are such that 
some of the mixing angles reach their upper experimental bounds. 
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4 Experimental bounds on sterile neutrino mixings 



The direct experimental searches for neutral leptons had been performed by a number of 
collaborations 28-4o|| (see e.g. 10, 13] for review of various constraints). The negative results 
of the searches are converted into the upper bound on $ a dp for different flavours. If neutrino 
oscillations are mediated by these sterile neutrinos, these bounds can be translated into the 
upper bounds on parameter \z\ and lower bounds on sterile neutrino lifetime. 

Below, we take a closer look at two main types of experiments ( "peak searches" and "beam- 
dump experiments" )0 and describe reinterpretation of these bounds in the case, when sterile 
neutrinos with MeV-GeV masses are also responsible for neutrino oscillations. 



4.1 Peak searches 



In "peak search" experiments [43| 
to charged lepton (e ± or /i ± ) and neutrino 
decay the limit on d 2 for masses in the range 60 MeV 



(see e.g. 



one considers the two-body decay of charged ti or K mesons 
10] for discussion). In case of the pion 
< M s < 130 MeV is provided by 
the searches for the secondary positron peak in the decay 7r + — > e + N to the massive sterile 
neutrino N as compared to the primary peak coming from the 7r + — > e + v e decay. Recent 
analysis of 39j puts this limit at d 2 < 10~ 8 in the mass range 60 — 129 MeV, for earlier results 



sec 



\m 3i 



30l . l3l| provided the 



In the smaller mass region (4 MeV < M s < 60 MeV) Refs. 
bound based on the change of the number of events in the primary positron peak located at 
energies M n /2 . Similar bounds were obtained for the same mixing angle in studies of kaon 



35 



decays [34| and for the in the decays of both pions [36j-[38| and kaons [34 , 

The lower bound on the sterile neutrino lifetime r s in the model (jl]), based on the peak 
search data and neutrino oscillations is shown in Fig. [1] by dot-dashed green lines. The 
parameters of neutrino mixing matrix are allowed to vary within their 3a limits (to minimize 
t s , while still keeping the values of all mixing angles compatible with the bounds from direct 
experimental searches). 



4.2 Beam-dump experiments and neutral currents contribution 



The second kind of experiments {"beam- dump experiments") [29j, |32j, |33j aims to create sterile 
neutrinos in decays of mesons and then searches for their decays into pairs of charged particles. 
Notice, that the expected signal in this second case is proportional to d\ or (and not to 
fl^ as in the case of peak searches, discussed in the Section H~T]) . We will demonstrate below 
that that in the models like (JH) (and in particular in the i/MSM) the results of some beam- 
dump experiments should be reinterpreted and will provide stronger bounds than discussed 
in previous works 10, [H, [TtJ (see also (i^). 



4.3 Reinterpretation of the PS191 and CHARM experiments 



The experiment PS191 at CERN was a "beam-dump" type of experiment described above [28 



29J. In searches for sterile neutrinos lighter than the pion M s < M^, the pair of charged par- 



8 Notice, that the neutrinoless double-beta decay (0vf3(3) does not provide significant restrictions on the 
parameters of the sterile neutrinos in the type-I see-saw models (contrary to the case discussed in e.g. Jioj), 
see discussion in 17, 41 1. In particular, this is the case in the i/MSM [42j . 
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Figure 1: The lower bounds on the lifetime of sterile neutrinos, responsible for the mixings 
between active neutrinos of different flavours in the see-saw models (HI). The bounds are 



based on the combination of negative results of direct experimental searches |28l . l29l . l3ll . 134- 
39 1 with the neutrino oscillation data (jjj]. The neutrino oscillation parameters are allowed to 
vary within their 3<r confidence intervals to minimize the lifetime. The solid black curve is 
based on our reinterpretation of PS191 data only, that takes into account charged and neutral 
current contributions (see Sec. 14. 3p . The interpretation of the PS191 experiment, taking into 
account only CC interactions (used e.g. in the previous works Jl3Ul7lP is shown in magenta 
dashed line. The bound from peak searches experiments only [3l|, I34j-l39| is plotted in green 
dot-dashed line. 



tides that were searched for in the neutrino decay comprised mostly of electron and positron: 



N 

e + e~v n 



(34) 



where N is a sterile neutrino with the mass M s . The first reaction in the chain is solely due 
to the charged- current (CC) interaction, and its rate is proportional to the ?9 2 . 

If sterile neutrinos interact through both charged and neutral currents (CC+NC) as it is 
the case in the models with the see-saw Lagrangian (j3J), any of three active-neutrino flavours 
may appear in the decay of N in (13~4"1) . The decay widths are 431 ]: 



T(iV -> eV 



967r 3 = 



with the following definition 



1 + 4 sin 2 9 W + 8 sin 4 9 W 



1 - 4 sin 2 6 W + 8 sin 4 6 



w 



(35) 



(36) 



and 9w is the Weinberg's angle so that sin 2 9w ~ 0.231 and c e ~ 0.59, c M ( T ) ~ 0.13. Therefore, 
the total number of events inside the detector that registers electron-positron pairs would be 
proportional to the combination of mixing angles $ 2 x (^ c a^a)- 



However, the model employed in the interpretation of the PS191 experiment [28|, |29 



was 



different, as has already been pointed in [44j. In the original analysis it was assumed that 



sterile neutrino interact only via charged currents, but not through neutral currents. In our 
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Figure 2: Comparison with the previous bounds on sterile neutrino lifetime in the z/MSM [17 



The solid purple curves represent the results of the present work, obtained by the combination 
of peak searches experiments 3jJ, I34M39I] together with the reanalysis of PS 191, that takes into 
account neutral currents (a union of black and green bounds from Fig. [1]). The red dashed 
curve is based on the combination of the same peak searches with the original interpretation of 
PS191 (i.e., with charged currents only). The blue dot-dashed line is taken from 171 ] . Notice, 
that the results of 17| were multiplied by a factor 2 to account for the Majorana nature of the 
particles (see discussion in Sec. I4.4j) . that was missing therein. The difference between the red 
and blue lines in the case of normal hierarchy is explained by wider 3a intervals for neutrino 
oscillation data, used in 



17] , compared to our work. 



language it means that 



he 



was used instead of the values flMM_l As was noticed 



above, the probability of meson decay into sterile neutrino does not alter if we exclude the 
neutral-current interaction, and therefore the total number of events with the electron-positron 
pair would be proportional to $1 x ■Q' 1 . 



Therefore if we denote the bounds listed in 
z/MSM takes form 



28. 29 




< ?9 4 

— e(errp) " 



then the bound for the 



(37) 



Similar bounds can be extracted from the reanalysis of meson decays into muon and sterile 
neutrino, that leads to replacement e — > /j, in (1371) . As a result, the reinterpretation of the 
results of the PS 191 experiment in combination with neutrino oscillation data produces up 
to an order of magnitude stronger bounds on lifetime than in the previous works (see Figs. [1] 
and E]). 

Similarly, the CHARM experiment [33] provided bounds on the mixing angles of sterile 
neutrinos in the mass range 0.5 GeV < M s < 2GeV. In the original analysis NC contributions 
were neglected. Therefore, to apply the results of this experiment to the case of the z/MSM, we 
reanalyzed the data as described above. In Fig. [3] we compare lifetime bounds coming from 
the CHARM experiment solely for CC and CC+NC interactions of sterile neutrinos. The 
difference in this case is about a factor of 20 



Model described in 28|, 2jJ contains only one Dirac neutrino, while in the ^MSM we have two Majorana 



fcrmions. Therefore actually c e = 1/2 in the original model. For details see Sec. 14.41 

10 In the case of the PS191 experiment, when using CC only for masses below the mass of pion suppression 
of the $g mixing angle due to neutrino oscillations meant that instead of -d\ bounds the lifetime is defined by 
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Figure 3: Comparison of the bounds on sterile neutrino lifetime (in the model 01])) based on 
the results of the CHARM experiments 33j solely (combined with the neutrino oscillation 
data). The orange (upper) curves correspond to the model with charged and neutral cur- 
rent interactions of sterile neutrinos, the brown (lower) - to the model with charged current 
interactions only. For details, see Sec. 14.31 



4.4 A note on Majorana vs Dirac neutrinos 

For completeness we briefly discuss the difference in interpreting experimental results for 
Majorana vs. Dirac sterile neutrinos. Similar discussion can be found e.g. in [ljjj]. When 
interpreting the experimental results one should take into account that in present work we 
consider two Majorana sterile neutrinos, while the experimental papers often phrase their 
bounds in terms of the mixing with a single Dirac neutrino, that we will denote U 2 . In the 
z/MSM twice more sterile neutrinos are produced per single reaction (because there are two 
sterile species - N 2 and N3), and, owing to their Majorana nature, each sterile neutrino decays 
twice faster (additional charge-conjugated decay modes are present). Notice, that the mass 
splitting between between two sterile states N 2 , N 3 is small \M 2 — M 3 \ <C \{M 2 + M 3 ) = M s 
and once born, the states oscillate fast into each other. Averaging over many oscillations can 
be accounted for by an extra factor | in the number of N 2 and N 3 species. Therefore, for 
beam-dump experiments one gets the same number of the detector events involving one Dirac 
sterile neutrino as one gets in the z/MSM if ("i^ + ^aa) 2 = Uts That is, one should identify 2 r d 2 a 
with the measured U 2 (recall fTTTj) that = |(^q 2 + ^03) )• I n the case of peak searches, the 
bound U 2 should be interpreted in the i/MSM as 2 + 3 <U 2 , as production of any state 
N 2 or iV^ contributes to the number of events in the secondary peak, i.e. again 2d 2 a should be 



identified with U 2 . Notice, that this factor 2 is missing in 17 . 



5 Results 

In this Section we summarize our results: the upper bound on the (combination of) mixing 
angles of sterile and active neutrinos in the see-saw models (j2J) in the range 10 MeV - 2 GeV 

the (much weaker) bounds. That led to the significant relaxation of the lower bound on the lifetime. If 
NC were taken into account, this was not possible anymore and therefore the lower bound on sterile neutrino 
lifetime became stronger by as much as the order on magnitude (black vs. magenta curve on the left panel in 
Fig. [1] In case of the CHARM experiment, both d 2 and d 2 ^ are strongly constrained and switching from one 
constraint to another makes (numerically) much smaller difference. 
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Figure 4: Direct accelerator bounds on the combination of active-sterile neutrino mixing 



angles, resulting from the reanalysis of the PS191 experiment |28|, l29j, taking into account 



decays of sterile neutrino through both charged and neutral currents and their Majorana 



nature. The shaded region is excluded. The case, analyzed in the original works [28|, |29 
(decay of sterile neutrino through the charged current only) corresponds to the choice c e = 1, 
c„ = c-r = 0, for details, see Sec. 14.31 We plot the bounds for two Majorana neutrinos (as in 



Fig. while in the original works [28l . |29| a single Dirac neutrino was analyzed. 



and the lower bound on sterile neutrino lifetime, obtained in combination of these bounds 
with constraints, coming from neutrino oscillation data. 



5.1 Bounds on the mixing angles of sterile neutrinos 

For the models (jlj) (two Majorana sterile neutrinos, interacting through both charged and 
neutral interactions), the compilation of constraints on various combinations of active-sterile 
mixing angles^ , i?^, # e Vz^ c a^, vS Ca ^a) that we used in this work are plotted in 
Figs. H and I51R 



5.2 The lower bound on the lifetime of sterile neutrinos 

The result of the Sections I3.2H3.31 combined with these experimental bounds can be trans- 
lated into the lower limits on the lifetime of sterile neutrinos. These results are presented in 
Figs. 6 on page 16 One sees that for the masses below that of pion the lifetime bound can be 
relaxed by as much as an order of magnitude (see the upper right panel in Fig. [H]). 



6 Conclusion 

In this work we have investigated experimental restrictions on the parameters of the see-saw 
Lagrangian in the case when two sterile neutrinos with the masses between ~ 10 MeV and 

11 Notice that in the published results of the PS191 experiment [29j bounds are given up to M s — 400 MeV. 
We extend these bounds up to 450 MeV, using the PhD Thesis of J.-M. Levy [iBl ]. 
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Figure 5: Direct accelerator bounds on the mixing angles. Left panel: -d 2 bounds, taken from 
3ll (region I), Q (region II), (reg ion III), (jjl-liill (region IV) and jiH (region V). Right 



panel: bounds, taken from |36H38| (region I), |34j (region II), [35j (region III), |29j] (region 
IV) and [32J (region V). The shaded regions are ruled out by the experimental findings. Dashed 
curve indicates mixing angle bounds given by original interpretation of PS191 experiment, but 
we do not use it to derive our final results, as explained in Sec. 14.31 The bounds are shown 
for the Majorana neutrino and are therefore two times stronger (see Section |4~4"|) . while in the 



original works [28|, |29[ a single Dirac neutrino has been considered. 



2 GeV, responsible for neutrino oscillations. Combined with the results of the direct experi- 
mental searches, the neutrino oscillation data provide stringent lower bounds on their lifetime, 
t s and allows to determine both maximum and minimum values of the mixing angles d 2 a . 



It was argued in 



171 ] that peak searches in the pion decay "might miss sterile neutrinos 
since the branching ratios can be extremely small due to the cancellation in d 2 , but the 
search of K + — > fi + + N can cover the whole allowed region" . As we have demonstrated 
in this work, the full suppression of d\ does not occur if the most recent data on neutrino 



oscillation experiments is taken into account. Moreover, contrary to the [l7j we conclude that 



the complete suppression of the mixing angle ft 2 may occur in the case of inverted hierarchy 
(Section 13. 3p . Our results imply that the peak searches in kaon decay may miss the sterile 
neutrino as well. 



We have reinterpreted the results of the PS191 experiment 28, 29], following 44] by taking 



into account not only charged-, but also neutral-current interactions (as both of these are 
present in the Type I see-saw Lagrangian). Our results demonstrated that in the low mass 
region (below the mass of the pion) the beam-dump experiments ($ 4 experiments) provide 
stronger restrictions than the peak search experiments (d 2 experiments) in case of normal 
hierarchy. In inverted hierarchy the reanalysis of the PS 191 experiment turns out to be very 
important as well. Also we point out that in the original analysis of the CHARM experiment 
32| neutral-current contributions were neglected as well. We have reinterpreted the results of 
this experiment in a similar way to PS191. The final results are presented in Figs. E 

Future experiments (for example, NA62 in CERN 1^1. LBNE experiment in FNAL^l or the 
LHCb experiment) have a great potential of discovering light neutral leptons of the i/MSM or 
significantly improving the bounds on their parameters (see e.g. 46[). 

^http : //na62 . web . cern . ch/M62| 
13 http : //lbne . f nal . gov 
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140 MeV 

Figure 6: The resulting lower bounds on sterile neutrino lifetime function of their 

mass, obtained by requiring that two Majorana sterile neutrinos are responsible for neutrino 
oscillations and their parameters do not contradict the negative results of direct experimental 
searches. In all figures the upper curve comes from using of the best fit neutrino oscillation 
parameters, the middle one - from their variation within the 3a limits, and the lower one does 
not take into account neutrino oscillation data and puts all three mixing angles equal to their 
direct experimental values. 
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Finally, we notice that sterile neutrinos with the masses in MeV-GeV re gion and small 
mixing angles can play significant role in the early Universe (for review see HI). Because 
of their small mixing angles these particles can be out-of-thermal equilibrium in the e arly 
Universe. In particular they can lead to the successful baryogenesis scenario [ill 
generate large lepton asymmetry at temperatures below the sphaleron freeze-out 
Fig. |7] we superimpose the bounds on \z\, coming from the direct experimental searches on 
the region of parameters (\z\, M s ) in which the successful bar yog enesis is possible (the region 




inside the black contours marked "BAU" based on the Ref. 24Q. Additionally, these sterile 
neutrinos can generate large lepton asymmetry at T < few GeV 23| that can affect production 
of dark matter sterile neutrinos 4fa, |49|. Their lifetime can be long enough to be present in 
primordial plasma at ~ MeV temperatures, so they can influence the primordial synthesis 
of light elements [12( and produce additional entropy and energy during their decays at later 
times 



50 



Sterile neutrinos with the mass of about 200 MeV and mixing angles ~ 10 — 10 
can affect the physics of supernova explosions Finally, our constraints may affect the low- 
reheating temperature cosmological scenario (see e.g. 52l. l53j|). 
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Figure 7: The region of successful baryogenesis in the z/MSM compared with the experimental 
lower bounds on the parameter \z\. The vales of M s and lying inside the black solid lines 
lead to the production of the observable baryon asymmetry (from fiih. The magenta dashed 
line marks is the lower bound on \z\~ 2 (parameter, called e in |23l. |24j) such that for smaller 
values at least one of the mixing angles d 2 a is in contradiction with direct experimental searches 
(under any assumptions about the values of unknown parameters of the PMNS matrix). 
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A Sterile neutrino lifetime 



For any phenomenologically interesting mass of the sterile neutrino it can decay to three 
neutrinos N — > vvv and it results in contribution 



vvv vvv vvv 



(38) 



to the total sum in Eq. (J2J) over decay products X. 

If M s > 2M e rj 1.0 MeV, then a new decay channel appears: N — > e + e~v. It corresponds 

to 



{- ±sin 2 9 W + 2 sin 4 9 W 



+ 12St (S e 4 - 1) In 
+ 2 sin 2 6 W (2 sin 2 9 W ±1 
+ 6S 4 e (l - 2S 2 e + 2St) In 



(1 - lASl - 2St - \2St) y/l - ASl 
l-3S 2 -(l-S 2 )v/r^4Sr 



S 2 e {2 + 105 2 - 12St)^l-ASl 
1-3S 2 -(1-S 2 )y^4gf 

52(i + yr^4f) 



where Sx = Mx/M s , 6w is the Weinberg angle, and sign "+" corresponds to a = e while "-" 
in the other case. For < \ we should add up to a sum terms that can be obtained 
from B° eu if one changes S e —> and if sign "+" corresponds to a — [A, "-" - to a = e, r. 
Other 5-coefficients are 



B$ = l-85j + 8Sj-Sj-24S>S' /il B^ = 0, (S, + S e <l), 

B« u = 6vr 2 ^(l - Slf, (S w < 1, M ff « 140 MeV, / ff « 130 MeV) , 
f 2 

B w = ^J^t 1 ~ ^) 2 > <^ M v~ 550 MeV, /„ « 155 MeV) , 



(39) 
(40) 
(41) 



Bl e = l^ 2 ^V 2 d {il-Slf-Sl{l + Sl)^ y/(l - (S w - 5 B )2)(1 - (S w + S e n (S e +S n < 1) , 

(42) 
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= 12tt 2 ^K 2 ,((1-^) 2 -^(1+^)) ^(l-^-S.mi-^ + S.n (S.+S^ < 1) . 

(43) 

In the last expression quantity V u d ~ 0.97 is the element of the Cabibbo-Kobayashi-Maskawa 
(CKM) quark matrix. Expression for B Ke , that is non-zero for mass range S e + Sk < 1 
(M K fa 495 MeV), can be derived from B e Ke by simultaneous change S n —> Sk, V u d Vus ~ 
0.23, f n — > fx ~ 160 MeV. Moreover, if we change in the resulting expression S e by S^, we get 
B Kii (for S^ + Sk < 1)- What concerns B-coefficients with a = r and X = (ne, iifi, Ke, Kfi), 
they are all zero. For masses of sterile neutrino, that are larger than the mass of the p-meson 
M p pa 780 MeV, other Bx appear. We do not list list them in the present paper, however they 
can be read from [13| . 

B PMNS parametrization 

In Section [3] we use the standard definition of the PMNS matrix (cf. e.g. [lj 

Sl3 \ / Cl2 S12 

V = | c 23 s 23 I I I [ -s 12 c 12 

1 

C12C13 C13S12 S13 

— C2ZS\2^ — C12S13S23 Ci2C23,e 1 ^ — Sl2Sl3 s 23 C13S23 (44) 
S23 s 12 e '^ ~ C12C23S13 — CyiSlZ^ ~ C23 s 12 s 13 C13C23 

where Cij = cos^-, and Sjj = sin^-. The active neutrino mixing matrix is given by expres- 
sion (J5]). 





C Ratio of sterile neutrino mixing angles for \z\ ~ 1 

As the expressions (1121) and f)20p show, the mixing angles have two terms: one is proportional 
to \z\ 2 and another to \z\~ 2 (recall that \z\ > 1). It was shown in Sec. 13.31 that for inverted 
hierarchy the |z| 2 -term can be zero for ft 2 and d 2 , while the \z\~ 2 term in general stays finite. 

For a given value of |z|, the |z| -2 -term is bounded from above. According to ffT2"]) and ([201) . 
its maximum is realized simultaneously with the maximal value of 

L N a H = \V a , + ^e ^ \[^V a2 \ 2 (45) 

V m 3 

in the normal hierarchy, and 

Li H = \V al + te^\^V a2 \ 2 (46) 
V mi 

in the inverted hierarchy. 

Analysis, similar to that of the Sections I3.2H3.3I shows that 

L N e H < 0.2, L" H < 1.1, L N T H < 1.1, L[ H < 1.96, L[ H < 1.3, L™ < 1.3 . (47) 
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These bounds allow to estimate the contribution of the \z\ 2 -terms to the whole sum of the 
squared mixing angles 



The ratio of (^ - (jigjl to (TJED gives 



T NH -i ri\ff i 
ffW = t RlH = L n f (AO) 

( 1+ mi)U|4 + 1 ' (! + 21) UI4 + ! ■ ^ 

V m 3 /II 1 v mi /II 

For z ~ 1 it can become of order unity. However, we restrict ourselves to the sufficiently large 
values z > 10, that are consistent with the upper bound, indicated by the experiments (see 
FigED 

i?f"<2xl(T 5 , <f<10" 4 , ^<2xl0- 4 , ^<10" 4 . (50) 

Comparison these results with the lower bounds (Table [2]) we see that z~ x terms are unimpor- 
tant for the for all mixing angles in NH and $ e in IH. What concerns the remaining angles $ M 
and d T in IH, they can be substantially modified by account of 2 _1 -terms, but anyway each of 
them can become small enough, compared to the other angles, as explained in next section. 
As a corollary, analysis and results of Sees. 13.2113.31 do not change significantly for large enough 
values of z. 
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